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Abstract 

Unitarity cannot be perserved order by order in ordinary perturbation theory 
because the constraint UU'^ = 1 is nonhnear. However, the corresponding con- 
straint for K = lnU, being K = —K'^, is hnear so it can be maintained in ev- 
ery order in a perturbative expansion of K. The perturbative expansion of K 
may be considered as a non-abehan generalization of the hnked-cluster expan- 
sion in probability theory and in statistical mechanics, and possesses similar 
advantages resulting from separating the short-range correlations from long- 
range effects. This point is illustrated in two QCD examples, in which delicate 

X" 

' cancellations encountered in summing Feynman diagrams of are avoided when 

they are calculated via the perturbative expansion of K. Applications to other 
problems are briefly discussed. 
Probability conservation is not maintained order by order in ordinary perturbation the- 
ory. This happens because the unitarity relation UU^ = 1 (for the time-evolution operator 
U) is nonlinear, whereas the constraints for the other exact conservation laws, such as energy, 
momentum, and charge, are linear. Order-by-order probability conservation can be restored 
if we expand instead K = lnU = J2n>i ^n/n, for then the unitarity constraint becomes the 
linear constraint K = —K'^. As long as every Kn is kept to be anti-hermitean, U = exp(i^) 
will be unitary no matter where the i^-expansion is truncated. We shall henceforth refer to 
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the perturbation theory in which K = \nU is expanded as the unitary perturbation theory. 
Mathematically, such an expansion is known as the Magnus expansion 

Besides unitarity, there is another reason to deal with the perturbative expansion Kn 
of K = InU, rather than U = 1 + J2n>iUn, because the former is the generalization of 
the linked-cluster expansion in probability theory and in statistical mechanics [Q, with the 
associated advantages in separating short-range correlations from the long-range effects. It 
is also analogous to the expansion of the effective potential in quantum field theory. The 
added complication here is that we are dealing with non-commutative operators rather than 
c- number functions. In QCD, short-range correlations refer to correlations in time and in 
color. This point will be explained later with two QCD examples. 

Unitary perturbation theory can also be regarded as an extension of the familiar eikonal 
approximation. For elastic scattering of a high-energy particle from a static potential, the 
scattering amplitude is diagonal in energy and the impact parameter. The scattering am- 
plitude {U) to all orders can be summed up to the form {exp{Ki)) = exp{{Ki)), where 
(Ki) is equal to the first Born approximation (Ui), and (• ■ ■) is the matrix element in 
impact-parameter space This eikonal approximation is also valid for electron-electron 

scattering when fermion loops are neglected |Q. For high-energy scattering of a quark 
from a color potential, amplitudes of all orders can again be summed up into the form 
{exp{K)) = exp{{K)), but now (K) is a color matrix, with (Kn) proportional to the nth 
nested commutator of the color matrices ta 0- In particular, (Kn) 7^ even for n > 1. 
Therefore unitary perturbation theory can be thought of as a generalization of the non- 
abelian eikonal approximation. 

The relation between f/„ and Kn can be obtained by expanding the exponential 
exp{J2 Kn/n) to the nth order. One obtains 

Un / ~r7TI~L. Kui-^KjYi^ ' ' ' KfYi,.! unk (1) 

K\ n,-i 

(m) •iij=i * 

where the sum is taken over all partition (m) = (mim2 ■ ■ ■mk) of the integer n. If U is the 
time-evolution operator between times T' and T, ordinary perturbation theory tells us 
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Un = {-iT I d^t Hj{h) ■ ■ ■ Hj{Q, (2) 

where TZn is the hyper-triangular integration region {T > ti > • ■ - t^ > T'}. Together with 
(|^), this may be used to derive a dynamical expression for Kn. In particular, 



Ki = Ui = -z f dt Hi(t), 

JT' 

K, = 2U2 - Ul = {-if r dh r dh [Hj{t,),Hj{t,)]. U12 (3) 

JT' JT' 

The difference between ordinary and unitary perturbation theories can be illustrated by 
a simple spin-^ example 0, with interaction —\(y ■B{t), and an external magnetic field 
B{t) = {Bj_cosujt, —B±smujt, Bo). The longitudinal field Bq causes a Zeeman splitting of 
magnitude Bo = Tiujq] the transverse field B± induces a transition between the lower and 
the upper states. Assuming the system to occupy the lower state initially, the probability 
to be in the upper state is plotted in Fig. 1, as a function of the scaled frequency variable 
Vt = h{uj — ujq)/Bj_ at the scaled time r = B±t/h = tt. The solid, deahed, dash-dot, and 
dotted curves are respectively the exact solution, the first Born approximation Pi = (f/i), the 
first unitary approximation UPi = {exp{Ki)) , and the second unitary approximation UP2 = 
(exp(i^i + K2/2)). The notation (■ ■ ■) represents the matrix element between the upper and 
the lower states. We see from Fig. 1 that the probability given by the Born approximation is 
larger than 1 for small violating unitarity, but both unitary approximations stay within 
the unitarity bound. As a function of r at the resonance frequency = 0, the exact solution 
as well as the two unitary approximations are periodic, with a period 2tt oscillating between 
and 1, whereas the Born approximation grows monotonically like r^, making it a worse 
and worse approximation at large r. 

For QCD, the unitary perturbation theory can be used to classify unitary parton-parton 



elastic amplitudes and to produce unitary model for Pomeron amplitudes [|TT|. 

Besides unitarity, it is often more profitable to compute Un via Kn, because Kn contains 
short-time and color correlations via the appearance of the commutator of the -ff/'s, instead 
of their product in Un- I will illustrate this remark with two second-order QCD examples. 



Higher-order situtations will be discussed later. 

The first example concerns high-energy electron-electron and quark-quark elastic scat- 
tering near the forward direction, with Mandelstam variables s and t. The only difference 
between the two cases is the presence of a color matrix ta at the vertices of QCD dia- 
grams, Fig. 2(a),(b),(c). The formulas for quark-quark amplitudes |g] are equally applicable 
to electron-electron amplitudes if we replace ta by 1. All amplitudes have a common fac- 
tor is/2m^ which will be factored out, leaving the rest to be (U). Fig. 2(a) gives the 
first Born amplitude (Ui) = {ig'^ /A'^)Gi, where A is the (transverse) momentum trans- 
fer with t = — A^, and Gi = ta x ta is the color factor of Fig. 2a. The amplitude for 
Fig. 2(b) is (-i^V27r)/2ln(se^^0G2, and that for Fig. 2(c) is {igy27r)l2\ns{G2 + GiNj2), 
where h = J d^k^/[{2nfkl{A - k^f]. The color factor G2 = Gi x Gi is that of 
Fig. 2(b), with Nc being the number of colors in the theory. This number is to be set 
equal to for QED. The second Born amplitude, being the sum of Figs. 2(b) and 2(c), is 
{U2) = (— 'i5'^/27r)/2(— 7rzG2 ~ In sGiNc/2). Note the cancelation of the Ins factor in the 
term proportional to 6*2- The Gi term consists of a color-octet exchange. It is not negligible 
compared to the first Born amplitude (Ui) in the leading-log approximation where (7^ ^ 1 
and In s = 0(1). In fact it is the beginning of the reggeized gluon contribution. The color 
factor G2 = Gi X Gi consists of two color-octet exchanges, which can be decomposed into 
a color-octet and a color-singlet. Its coefficient in the impact-parameter representation is 
equal to the coefficient of Gi x Gi in \ {Ul), a factor of Ins down from the coefficient of the 
Gi term. Its color-octet contribution is therefore negligible in the leading-log approxima- 
tion, but its color-singlet component is leading and constitute the beginning of the Pomeron 
amplitude. For QED, the factor in front of Gi is zero, and G2 is replaced by 1. Hence 
(f/2) = in the impact-parameter representation. 

To summarize, a delicate cancelation of Ins occurs between Figs. 2(b) and 2(c) in QED, 
and in the color-singlet channel of QCD. This is reminiscent of the delicate cancellation of 
the volume factor in the grand partition function in statistical mechanics unless the linked- 
cluster expansion is used 0. As a result, the k color-octet exchange amplitude in QCD 
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is proportional to g'^^ {g'^ In s)'^~'' {k = 1,2), and it is this dependence that leads to the 
reggeon amplitudes in QCD. We shall now see that if we calculate {U2) via (Kf) and {K2), 
then delicate cancelation is not needed, and the energy dependence needed for the reggeon 
structure becomes immediate. 

Since {Ki) = (Ui), it is given simply by Fig. 2(a). For {K2/2), it can be computed 
from Fig. 2(c), provided the product of the color matrices tatt in the upper quark line is 
replaced by their commutator [ta, U] p, thus making the color factor purely Gi. Otherwise, 
the result is identical to that of Fig. 2(c). In the impact-parameter representation where 
Ki are diagonal, we obtain in this way the same result as before: {K2/2) is given by the 
Gi coefficient of Fig. 2(c), or that of (f/2), and {Ki)"^ /2 is given by the G2 = Gi x Gi 
coefficient of (f/2)- The absence of the Ins factor in {Ki)'^/2 is not due to cancelation; it 
is simply a consequence that {Ki) is independent of s. This is the same as in statistical 
mechanics when linked cluster expansion is used 0. The reggeized gluon contribution is 
now isolated in {K2/2) alone, without having to combine two Feynman diagrams to obtain 
it. Its octet origin can be traced back to the color commutator in the upper quark line, or 
the commutator structure of K2 shown in (^. This gives {K2) a physical meaning as (part 
of) the one-reggeon amplitude. 

The second example is pion-nucleon elastic scattering in the large- A^^c limit. The tree 
diagrams are shown in Fig. (d),(e), with pion-quark interaction given by a Hamiltonian of 
the form Hi ~ g^^F-vrg, where F is a matrix that contains the isospin and spin information. 
The effective coupling at each vertex is proportional to {Hj)/^/Wc, where (■ ■ ■) represents 
the nucleon matrix element and l/y/TT^ is the normalization factor needed for each external 
pion. Since a color-singlet nucleon contains Nc quarks, (Hi) ~ Nc, so the effective coupling 
at each vertex is of order y/N^, making each tree amplitude ~ N^. However, the term 
proportional to A^^^ is canceled in the sum of the two diagrams, leaving a total amplitude 
(f/2) of order unity in the large-A^c limit. This can be understood ^ directly from (^ 
firstly because (Ki) = on account of energy- momentum conservation: an on-shell nucleon 
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cannot absorb or emit a massive pion and remains on-shell. Thus (f/2) = \{K2). Since 
the commutator of two one-body operators is again a one-body operator, hence ([-f/^i, ^^2]) ~ 
(g^[f i-TTi, f'2-T^-2\q) ~ N(.. Taking into account the normahzation factor 1/ \/Wc for two pions, 
we conclude that (f/2) = \{K2) ~ N^/Nc = 1, as needed. Like the situation in Example 1, 
again there is no need to cancel the Nc term explicitly, because the commutator structure 
of K2 already provides for it. Physically, the one-body nature of K2 tells us that the two 
pions must interact with the same quark in the nucleon. 

With these two examples, it is clear that the commutator structure of K2 is the key to the 
simplification. The question is whether Kn still possesses such commutator structures for 
higher n. The answer is 'y^s', though the detailed structure is increasingly more complicated. 
Let us define anti-hermitean operators C„ to have the simple nested commutator structure, 

Cn=Hrf d^t[H,,[H2,[---,[H^.uHn] ■■■]]], cn (4) 



where Hi = Hj{ti). Then it can be shown [jlO| that 



k ^ 

Un ^ ^ W CmiCfn2 ' ' ' Cmi^i xic (5) 

(m) j=l ^i=j 

where the sum is taken over all partitions (m) = (mim2 ■ ■ ■ ruk) of the number n. In 
particular, Ui = C, = K,, U2 = {Cj + C2) /2 = {Kf + K2) /2, f/3 = {Cf + C1C2) /6 + 
(C2C1 + C3) /3 = Kf/Q + {KiK2 + K2K1) /4 + i^'3/3. From these relation, or more generally 
by comparing (|l|) with (|^), one obtains the commutator structure of to be = Cn + Rn, 
with Ri = R2 = 0, i?3 = [C2,Ci]/4, and R^ = [C3,Ci]/3. In general, i?„ is given by 
commutators of Cm, with increasingly complicated coefficients as n increases. 

We can use these relations between f/„ and Cm to calculate quark-quark and electron- 
electron scatterings to higher orders, avoiding delicate cancelations, and obtaining directly 
the energy dependence necessary for reggeon structures in QCD, viz., an amplitude propor- 
tional to g^^{g'^\Yis)'^~^Gk when k color-octet objects are being exchange in the t-channel 
PJTT|]. They can also be used to obtain the correct amplitude |T2| ~ A^^^~"/^ for pi-nucleon 



inelastic scattering with n — 1 pions in the final state, without having to deal with delicate 
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cancelation of n — 1 powers of Nc encountered by summing Feynman tree diagrams . As 
mentioned at the beginnig, abelian and non-abelian eikonal approximations are special cases 
of these formulas. Infrared structure of QED can be obtained from these formulas, with Ki 
containing the Bloch-Nordsieck result, K2 giving rise to the Coulomb phase, and Kn = for 
n > 3 |TB|,|T^. The Landau-Pomeranchuk-Migdal effect \T5\ can also be cast into this form 
|T6| with Kn = for n > 3. Other applications may require knowing the the diagrammatic 
rules for calculating Kn- For QED such rules are known fl^. For tree diagrams in QCD 
involving an energetic particle such rules are also known [^, but the general case for QCD 
is still to be worked out. 

So far we have emphasized situations when U is the time-evolution operator with a 
known Hamiltonian. If the dynamics is unknown, we can still use Kn in (|l]), or C„ in (^, 
to parametrize the unitary dynamics. As long as Kn is anti-hermitean, it is clear that U 



is unitary. As long as Cn is anti-hermitean, it can also be shown |jT^ that U is unitary, 
though the proof is much more involved. A Wolfenstein-like parametrization of the CKM 
matrix to all orders can be deduced from these parametrizations. They can also be used to 
study the unitary matrix describing the overlap of the unperturbed and the perturb energy 
eigenfunctions [|T^ . 

In conclusion, we have given a number of examples to show that the unitary perturbation 
theory, in which K = Inf/ is expanded perturbatively, is in many ways better than the 
ordinary perturbation theory, where U is directly expanded. 

This research is supported in part by the Natural Sciences and Engineering Research 
Council of Canada, and the Fonds pour la formation de Chercheurs et I'Aide a la Recherche 
of Quebec. I am indebted to Saad Ali, Hong- Mo Chan, Marc Grisaru, Greg Mahlon, Sheung- 
Tsun Tsou, Wu-ki Tung, and Kenneth Young for stimulating discussions. 
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FIG. 1. Unitary perturbation theory vs ordinary perturbation theory in a two- level example. 



(a) 



(b) 



(c) 



(d) 



(e) 



FIG. 2. (a,b,c), high-energy quark-quark scattering diagrams to one loop order. (d,e), me- 
son-nucleon scattering at large Nc in the tree approximation . 
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